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The semiclassical Wigner treatment of Brown and Heller [J. Chem. Phys. 75, 186 (1981)] is applied 
to triatomic direct photodissociations with the aim of accurately predicting final state distributions 
at relatively low computational cost, and having available a powerful interpretative tool. For the 
first time, the treatment is full-dimensional. The proposed formulation closely parallels the quan- 
tum description as far as possible. An approximate version is proposed, which is still accurate while 
numerically much more efficient. In addition to be weighted by usual vibrational Wigner distribu- 
tions, final phase space states appear to be weighted by new rotational Wigner distributions. These 
densities have remarkable structures clearly showing that classical trajectories most contributing 
to rotational state j are those reaching the products with a rotational angular momentum close 
to + l)]^''^ (in n unit). The previous methods involve running trajectories from the reagent 
molecule onto the products. The alternative backward approach [L. Bonnet, J. Chem. Phys. 133, 
174108 (2010)], in which trajectories are run in the reverse direction, is shown to strongly improve 
the numerical efficiency of the most rigorous method in addition to be state-selective, and thus, 
ideally suited to the description of state-correlated distributions measured in velocity imaging ex- 
periments. The results obtained by means of the previous methods are compared with rigorous 
quantum results in the case of Guo's triatomic-like model of methyl iodide photodissociation [J. 
Chem. Phys. 96, 6629 (1992)] and an astonishing agreement is found. In comparison, the standard 
method of Goursaud et al. [J. Chem. Phys. 65, 5453 (1976)] is only semi-quantitative. 



I. INTRODUCTION 

Modern molecular beam and spectroscopic technics allow the measurement of quantum state distributions of pho- 
todissociation products with an amazing level of precision [iHS] . Accurate theoretical descriptions of the mechanics 
of molecular fragmentation 6J are thus needed to reproduce and rationalize these data, or predict them when ex- 
periments cannot be performed. Beyond their fundamental interest, photodissociation dynamics studies provide very 
useful data to specialists of planetary atmospheres [7j or interstellar clouds [8], and they form a branch of molecular 
physics which has continuously generated technological progress for more than a century. 

The goal of the present paper is to report some new insights into the semiclassical Wigner treatment [SI llOj of the 
dynamics of direct triatomic photodissociations (see also refs. [6| IlllfIS] for related works). In these very common 
processes, nascent products strongly repel each other after the optical excitation and are completely free in a few 
tenths of femtoseconds [5]. The practical interest of the semiclassical Wigner method, proposed about three decades 
ago by Brown and Heller [3], is that it was found to be more accurate than the earlier and more classical approach of 
Goursaud et al. O [H] , and leads to calculations much easier to perform than exact quantum calculations jTOHST] . As 
far as polyatomic molecules are concerned, the latter are usually prohibitive and the semiclassical Wigner treatment 
might thus be an interesting alternative. From the fundamental side, the interpretative power of this approach makes 
it a powerful tool for rationalizing the dynamics. 

The semiclassical Wigner method is based on the notion of Wigner distribution associated with a given quantum 
state [21 l32H39j . For a one-dimensional system of configuration coordinate x and conjugate momentum px in the 
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state l^*), the Wigner distribution is defined as 

p(x,p^)^^f dse^'P-'/'' (^\x + s) (x~ s\^). (1) 
ttH J 

This fascinating density was introduced by Wigner in 1932 [32 as a convenient tool for quantum mechanically 
correcting the Gibbs-Boltzmann theory of thermodynamic equilibrium in the limit of small temperatures. Integration 
of p{x,Px) over leads to Kxl^*)]^ while integration over x leads to p. In addition to that, calling respectively 

H and H the classical and quantum Hamiltonians of the system, integration of p{x,Px)H with respect to x and Px 
leads to One might thus be tempted to call p{x,px) the phase space distribution corresponding to l^*). 

Stricto sensu, however, one should not, for p{x,px) takes negative values in some parts of the phase space when 
{x\^) (or {px\^)) involves at least one node, and the uncertainty principle puts a shadow on the phase space concept. 
Nevertheless, one can formally use this non conventional density of probability as an usual one. 

The beauty of the semiclassical Wigner treatment is that it mixes the quantum and classical descriptions in a very 
natural way [21 [TU]. The underlying principle of the treatment is as follows. The population of a given product 
quantum state is shown to be proportional to the square modulous of the overlap between (i), the product state 
wave function, and (ii), the time-evolved wave packet resulting from the propagation of the initial reagent state 
(multiplied by the transition moment) over a period of time large enough for the whole wave packet to be located 
in the product channel [5J [5] . This expression is then transformed into an overlap between the Wigner distribution 
corresponding to the product state, and the Wigner distribution associated with the time-evolved wave packet [9]. 
Everything is rigorous up to now. What makes the Wigner treatment semiclassical is that the Wigner distribution 
of the time-evolved wave packet is obtained by propagating in time the Wigner density of the initial reagent state 
according to the laws of classical mechanics 9J, just as if this density was a solution of Liouville equation [40]. To 
recap, initial and final dynamical states (phase space points) are assigned statistical weights according to quantum 
mechanics through Wigner distributions while nuclear dynamics are ruled by classical mechanics. The corresponding 
mathematical developments are given in section [II A| In this approach, trajectories are run forward in time from the 
reagent molecule onto the separated products. We shall call it forward I. 

In the only two applications of the semiclassical Wigner method that seem to have been published to date [SI dO] j 
rotation motions were frozen. In the present work, however, we take into account the full-dimensionality of the 
triatomic system, making thereby the method applicable to processes taking place in realistic conditions. 

In addition to the usual vibrational Wigner distributions [HI UHl UH UHl [55Vl38j . the method involves rotational 
Wigner distributions which, to our knowledge, are introduced for the first time in the context of reaction dynamics. 
These quantities are respectively discussed in sections [11 B| and p^I C| Rotational Wigner distributions have remarkable 
structures clearly showing that the classical trajectories most contributing to rotational state j are those reaching the 
products with a rotational angular momentum close to + 1)]^^'^ (in fi unit). 

The study of a Franck-Condon model process in section |IID| allows to illustrate how the degrading effect, discussed 
by Gray and Truhlar [13] and Schinke [10], alters state-resolved cross sections en route to products. This suggests a 
slight formal modification of the forward I method in order to greatly improve its accuracy. 

Nevertheless, this method turns out to have a limited numerical efficiency. An approximation is thus proposed in 
section [II E| in order to strongly increase it. The resulting method will be called forward II. 

The recent backward approach |16H18j . in which trajectories are run from the products onto the reagent molecule. 



is shown in section III to be strictly equivalent to the forward I method while being numerically much more efficient. 
In addition to that, it is state- selective, and thus, ideally suited to the description of state-correlated distributions 
measured in velocity imaging experiments [iHS]. 



Since we shall also apply the standard method of Goursaud et al. [IT] , we briefly recall its main lines in section IV 
The results obtained by means of the two forward methods and the backward one are compared in section [V] with 
rigorous quantum results [221 130] in the case of the reduced dimensionality model of methyl iodide photodissociation 
of Guo (4TI, and very good - often quantitative - agreement is found. In comparison, the method of Goursaud et al. 
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[TT] is only semi-quantitative, at least for one of the two excited electronic states involved in the process, within which 
the system has a more quantum-like behavior. In particular, this approach does not allow to systematically reproduce 
inverted vibrational or rotational state populations, contrary to the semiclassical Wigner treatment. Section |VI| 
concludes. 



II. FORWARD SEMICLASSICAL WIGNER APPROACH 
A. Accurate formulation 

1. System 

Let us consider a triatomic molecule ABC optically excited by a photon of energy hv from its rovibronic ground 
state at energy Eq up to a given repulsive electronic excited state. The future products, say AB and C, strongly 
repel each other immediately after the photon absorption and are formed in a few tenths of femtoseconds. The 
dissociation is thus direct. The total energy of ABC \s E — E^) + hv. The optical excitation is supposed not to 
excite the rotation motion. Within the framework of this reasonable approximation, the total angular momentum J 
is kept at throughout the whole fragmentation process (see refs. [HI UHl US 1131 HI] for rigorous treatments of the 
possible rotational transitions implied by an optical excitation). R is the vector going from the center-of-mass G of 
AB to C and r is the one from A to B. Vg and Vg are the potential energies in the ground and excited electronic 
states, respectively. They both depend on the moduli R and r of R and r and the angle 6 between them. Far in the 
asymptotic channel, Ve reduces to the potential energy of the free AB diatom, denoted Ve{T). The zero of energy is 
defined as the minimum of Vg, or Ve{r), in the separated products, m and are, respectively, the reduced masses of 
AB, and C with respect to AB. d = d(R, r) is the transition dipole vector responsible for the electronic transition [S], 
and e is a unit vector in the direction of the polarization of the electric field of the photon. The final quantum state 
of AB is denoted (n,j), where n and j are the vibrational and rotational quantum numbers, respectively. 

2. Quantum partial cross section and product state distribution 

Within the electric dipole approximation, the (n, j) state-resolved, or partial, absorption cross section is propor- 
tional to |6] 



*B'|d.e|</.o) 



(2) 



l^'^-') is the state of inelastic scattering between AB and C at energy E, in the electronic excited state, and with 
outgoing free wave in channel {n,j). |0o) is the rovibrational ground state in the electronic ground state. Final state 
populations are deduced from Eq. ([2| according to 



A very common approximation consists in replacing in Eq. ([2| d.e by the modulous d of d, which amounts to suppose 
that d is parallel to e at the instant of the photon absorption. This approximation appears to be very satisfying 
provided than one is not interested in vector properties IMl UHl HTHi5] . d is either parallel or perpendicular to the 
triatomic plane and d does only depend on the configuration of ABC, i.e., d = d{R,r,9). 
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Setting 



3. Quantum partial cross section at time t 



Eq. (|2| reads 



Writing the identity operator as 



$o(R,r) =d(i?,r,0)4(R,r), 



/ 



dRdr ^'^^ *(R,r)$o(R-, r) 



(4) 



(5) 



(6) 



where the optical excitation defines the origin of time, and inserting this operator between the two states in Eq. ([s]) 
leads to 



J dRdr *^-'*(R,r)e'^*/'^l>t(R,r) 



(7) 



The phase factor e'-^*/^ arises from the fact that l^'^'') is a stationary state of the system in the electronic excited 
state. This factor is only written for clarity's sake, since its modulous is equal to one. $t(R, r) is the wave-packet 
obtained by propagating ^oCR-jT) in the electronic excited state during t. 
It is shown in Appendix A that Eq. Q can be rewritten as 



dRdrdO sine ■^"J*iR,r, 



(8) 



where ^'^■''(i?, r, 6*) = 8^ / ^ n Rr^/"J (R^r) and <S>t{R,r,e) = 8^/^TrRr^t{'R,r). ^^-^ (i?, r, 61) and $f(i?,r,6') are shown in 
the same appendix to satisfy the Schrodinger equations 



(9) 



and 



with 



ih 



dt 



T ^ — " ' — SITlu 

2^ai?2 2m dr"^ 2Isin0 80 89' 



(10) 



(11) 



/ is the reduced moment of inertia given by Eq. (A. 12) 



^t{R,r,9) is the wave-packet obtained by propagating ^o{R,i^,S) during t. Setting <j)o{R,r, 
we deduce from Eq. Q and Eq. (A.5) at i = 0, the identity 



= 8i/27ri?r(/)o(R,r), 
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^o{R,r,e)^diR,r,e)MR,r,e). 



(12) 



bo{R,r,6) is solution of 



Hgcl)o = [f + Vg]<l>o = Eo<l>o- 



(13) 



Far in the asymptotic channel, the general expression of VE*^^ {R, r, 9) is 



*^^(i?,',> 



-,1/2 



27r?i kn 



n'i' 



-,1/2 



27r?i fc„'-,' 



e-^^-'^'\n'{r)Yf.{e). (14) 



This expression assumes that the vibration and rotation motions of AB are uncoupled, a good approximation 
provided that too highly excited rovibrational states are not available, knj is defined by 



fc„j = ^[2/i(^-£;,y)]'/' 



(15) 



where Enj is the internal energy of AB in state {n,j). Xn(^) is the nth excited vibrational state, and Yj\6) is the 
jth spherical harmonic defined by 



2j + l 



An 



1/2 



Pj (cosk 



(16) 



Pj{x) is the jth Legendre polynomial. The S'-matrix element Snjn'j' is the probability amplitude to go from {n',j') 
to {n,j). The ^^-^ (i?, r, ^?)'s satisfy the usual orthogonality relations. 

For t tending to +oo, $t(i?, r, 0) entirely lies in the product channel and is moving outwards. The overlap in 
Eq. ([8| between the incoming part of '^^l^ {R,r,9) (see Eq. (14)) and $i(i?, r, 0) is thus zero. Hence, Eq. ([s]) becomes 



lim 



dRdrdO sinO 



-,1/2 



e-^'=-V(r)Y,"(0) $t(i?,r,0) 



(17) 



{Xn{f) and Yj{6) being real functions, they are equal to their complex conjugate). 



4- Quantum partial cross section m terms of Wigner densities 



It is shown in Appendix B that a strictly equivalent phase space integral expression of Eq. (17) is 



E'^J' = lim {2^nf [ dr ptiT) ptr{R,P) Pn{r,p) Pj{e,Pe) 



(18) 



with r = {R,r,9, P,p, Pe). pt{T) is the Wigner density related to ^t{R,r,9) through Eq. (B.3). ptr{R,P), Pn{r,p) 
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and pj(6,Pe), respectively called translational, vibrational and rotational Wigner distributions, are given by 



1 r ^ k] 



2fi 2/i 



e(P), 



(19) 



and 



Pnir,p) [ ds e'^P^/^ Xn{r + s) x„(r - s) 



pAe,Pg) ^ \ ( ds e^'^"/'^' sin{9 + s) ¥^{9 + s) sinlO - s) Y"{9 - s). 



(20) 



(21) 



The present expression of ptr{R,P) (Eq. (19)) appears to be different from the one in refs. [9l[T^, due to different 

and 



normalizations of translational states in Eq. ( 14). p„(r,p) and pj{9, Pg) are discussed in sections 



II B 



lie 



The ar- 



gument of Y,° being limited to the range [0,7r], the constraints on s in Eq. (21 ) are —9 <s<tt — 9 and 9 — ii<s<9,or 



max[—9, 9 ~ t:] < s < min[9, n — 9]. 



(22) 



5. Passage to the semiclassical description 

We now introduce in the previous rigorous quantum formulation the following classical ingredient: we assume that 
the propagation of Pr(r) from r = to r = i is correctly ensured by classical mechanics, i.e., Pri^) reasonably 
satisfies Liouville equation [ID] . As is well known, this assumption is only valid over a short period of time [TUl [131 155] 
and its application should thus be limited to direct photodissociations. We shall come back to this important issue 
later below. 

In the framework of the previous assumption, we have according to Liouville theorem |40] 



Pt{T)dT = po(ro)dro. 



(23) 



In this identity, T should be understood as the dynamical state of ABC reached at time t when starting from state Fq 
at time 0. For both clarity and simplicity's sake, the components of F will be called {Rt,rt,9t, Pt,Pt, Pet) in the rest 
of this section. In practice, they are determined by solving the Hamilton equations given in chapter 5 of Schinke's 
book (see Eqs. (5.3) and (5.4) where to, /i, V and 7 are respectively denoted here by fi, m, Vg and 9 ; these equations 
are clearly the classical analogs of Eqs. ([9| and (11) in the present work). Hamilton equations involve the classical 
Hamiltonian 



p2 

2p 



P 

2m 



Pi 
21 



Ve{R, r, 9) 



(24) 



needed in the followings. 



Eqs. (18) and (23 1 finally lead to the semiclassical Wigner expression of Brown and Heller fTU] 



lim {2nhf [ dTo po{To) PtriRt,Pt) Pn{rt,Pt) Pj{Ot,Pet) 
^+00 J 



(25) 
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except that here, the rotation motion of AB is taken into account. 

A comment on 0t is in order. As previously seen, Y^{9t) is only defined for 6t in the range [0,7r], and we shall see 
later below (section IIC) that this is also the case of pj{6t, Pet)- However, 9t is a (roughly) linear function of time 



in the asymptotic channel (see section II D ) and thus, it eventually gets out of the previous interval. Replacing 6t in 



Eq. (25) by 



TT int(9t/n) 



(26) 



if int{0t/n) is even, or 



Tr~[0t~7Tint{0t/Tr)] 



(27) 



if int{0t/n) is odd, keeps the ABC configuration unchanged and constraints to belong to the range [0,7r] 



Nevertheless, we shall see in section II D that Eq. ( 25 ) , despite its solid theoretical grounds and its elegance, cannot 



be applied as such. This is due to the degrading effect, consequence of the limitation of the validity of Eq. (23) to 



short times. Eq. (25 1 needs a slight formal modification, proposed after discussing to some extent the degrading effect. 



B. Vibrational Wigner densities 

Vibrational Wigner densities have been discussed elsewhere, in particular for the harmonic and Morse oscillators 
[51 I37[ I38| . Therefore, we only concentrate on those involved in the process studied later in this work in order 
to check the validity of the semiclassical Wigner method, i.e., the reduced dimensionality model of methyl iodide 
photodissociation of Guo [41] . 

In this model, the methyl radical vibration motion is reduced to its umbrella mode, treated as if this was the 
stretching mode of a pseudo diatomic molecule. The coordinate r is between the centers-of-mass of the three hydrogen 
atoms and the carbon atom, assuming the C^y symmetry is preserved throughout the whole process (see section [v] for 
more details). The potential energy i'e(f), represented in Fig. [l] together with its second order development, clearly 
appears to be strongly anharmonic with a significant contribution from a fourth order term. 

The first vibrational states have been calculated by means of the Truhlar-Numerov algorithm [33] and perspective 



views of the resulting Wigner distributions, estimated from Eq. (20) over a regular grid of 100 points, are displayed 
in Fig. [2] for the levels n = — 5. The r and p axis are directed towards the right and left, respectively, r belongs to 
the range [-1.2,1.2] and Pr to the range [-12,12], both in atomic units. Despite the strong anharmonicity of Ve{r), the 
shape similarity with the Wigner distributions of the harmonic oscillator is stricking [6] |38] . 

While for n — 0, the Wigner density is always positive, it takes negative values for n > 1 in the crater-like areas. 
Craters are delimited by cones on the edge of which small summits are visible, particularly for n > 3. Contrary to 
the harmonic case, contour levels do not exactly correspond to classical orbits, especially in the vicinity of the edge 
[45] . This is at the origin of the degrading effect [TOl [13] , as discussed further below. This effect will however appear 
to be much stronger in the case of the rotation motion. 



C. Rotational Wigner densities 

Perspective views of pj{0,Pg), estimated from Eq. (21 1 over a regular grid of 100 points, are shown in Fig. |3] for 
j — 3, 6 and 10. The and Pq axis are oriented towards the right and left, respectively. belongs to the range 
[0,7r], and Pe to the range [-15?i,15?i]. 

At first sight, pj{0,Pe) appears to involve j narrow wells along the symmetry axis of the distribution defined by 
Pq = 0, separated by j — 1 peaks (the second part of this statement is not true for j — Q). These wells are in fact 
negative peaks with magnitudes comparable with the ones of positive peaks. This is clearly seen in Fig. [4] where a 
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front view of pQ{9,Pe) along the ^-direction is displayed. In this complementary view where the previous peaks are 
aligned, positive ones are in the light while negative ones are in the shadow. We shall call central peaks this chain of 
alternatively positive and negative spikes. In addition, pj(9,Pg) for j involves two positive lateral ridges parallel 
to the 6'-axis (note that for j = 0, the two lateral ridges have merged with the central peak). As indicated in Fig. |4j 
the summits of these ridges, located along the line 9 — n/2, appear to be roughly defined by Pg « ±h[j{j + 1)]^^^ in 
the upper and lower half plane, respectively. This is actually true only for large j'a, as illustrated in Fig. [5] where the 
difference between (i) the exact value of Pg corresponding to the summit in the upper half plane and (ii) the quantum 
value h[j{j + 1)]^/^ is seen to decrease to in terms of j. 

For completeness, the semiclassical limit of the rotational Wigner density is considered in Appendix C. The interest 
of this limit is that it clearly explains the main topological features of pj{9, Pg) outlined above. Moreover, it provides 
analytical expressions which prove to be useful to make the link between the semiclassical Wigner treatment and the 
method of Goursaud et al. 11;. The study of this link will be published elsewhere |46j . 

Due to the sign alternation of central peaks, their contribution to E^-*, given by Eq. (25), is expected to be negligible 
as compared to the one of lateral ridges provided that the distribution of the points {9t,Pgt) overlaps several peaks. 
In the semiclassical limit of large j's where the central peaks become very sharp (see pio{9,Pg) in Fig. [s]), this will 
necessarily be the case. Consequently, one recovers the well known semiclassical rule stating that the values of Pg 
close to ±h[j{j + 1)]^^^ ~ ±h{j + 1/2) mostly contribute to the population of the j*'* rotational state. Conversely, 
the contribution of central peaks to for j = 1 — 3 cannot be excluded and the semiclassical rule fails (in the case 
where j = 0, the semiclassical rule works for the trivial reason that there is only one central massif). 

In the next two paragraphs where vibrational and rotational Wigner distributions are compared, we use the quantum 
number j for both densities (and for the vibrational eigenstate as well). The vibrational density pj{r,p) and its 
rotational analog pj{6,Pg) are then differentiated from their arguments only. 

The functions sin9 Yq{9) and sinO Y^{9) are proportional to sin6 and sin9cos9, respectively. Topologically, these 
two functions are not very different from the ground and first excited vibrational states of a diatom. This is the reason 
why they result in densities pa{9,Pg) and pi{9,Pg) bearing strong similarities with pQ{r,p) and pi{r,p), respectively. 
This is clearly seen when comparing Figs. [2] and [3j For the ground states, both densities are bell-shaped and for the 
first excited states, they both involve a crater with small summits lying on its edge (smaller for the vibrational density 
than for the rotational one). 

On the other hand, sind Yj'(9) is significantly different from Xj(^) for J ^ 2, in particular because of the sin9 
term which makes the oscillation amplitude of sin9 Yj'{9) decrease when going away from 7r/2 (see Appendix C) 
while the one of Xji''') tends to increase when going away from the equilibrium geometry (at least up to the classical 
turning points). Consequently, pj{9,Pg) appears to be very different from pj{r,p) for j > 2, as can be observed when 
comparing Figs. [2] and [3j 

For the free rotor, Pg is a constant of motion and classical orbits are defined by straight lines parallel to the 0-axis. 
It is thus clear from Figs. |3] and [lO] that the rotational Wigner density strongly varies along a classical orbit, a fact 
responsible for the degrading effect, as shown in the next section. Note that the amplitude of the previous variation 
is much stronger than for the vibration motion. On average, classical orbits may indeed be shown to be much closer 
to contour levels for the vibrational Wigner density than for the rotational one. 

The Wigner distribution of a rigid rotator has already been discussed in ref. [17]. However, the context is quite 
different, hence leading to a different mathematical definition of the distribution. In addition, it seems that the notion 
of rotational Wigner density has never been introduced in the reaction dynamics field to date. Lastly, we note the 
shape similarity between the present rotational Wigner distributions and those for a symmetric infinite square well 
potential (compare Fig. 1 in ref. [SB] and Figs. |3] and 10 in the present work). 
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D. Franck-Condon process 



1. System 



We now assume that Vg is isotropic, i.e., does not depend on 0. Moreover, we freeze the AB vibration motion for 
simphcity's sake. AB is thus a rigid rotor the length of which is denoted r,,. In addition, the reduced moment of 



inertia /, given by Eq. (A. 12 1, is supposed to reduce to mr^. This is a good approximation, since generally, R is 



already larger than Tg at time 0, and unless C is much lighter than both A and B, /i is also larger than m. Lastly, 
the transition dipole moment is kept at a constant value. 

2. Degrading effect 



Following the developments of section II A the j state-resolved absorption cross section is proportional to 



y — 



dRdO sine l-^ {R,e) <i>o{R, 



(28) 



(see Eq. ([s]) at time without the r coordinate). Due to the isotropy of Ve which makes the radial and angular 
motions uncoupled, the scattering state \I''^(i?, 9) can be written as 



*i,(i?,0)-C/B(i?)Y,°(0). 



(29) 



Asymptotically, 



Ue{R) 



1/2 



(30) 



with 



2/i E 



^^jU + 1) 

2TOr? 



1/2 



(31) 



r] is the phase shift. We shall suppose that T4 is sufficiently repulsive for the final translational energy to be much 
larger than the rotational energy h?' j{j + 1) / {2mr'^). In other words, the product energy E is mainly deposited into the 
translation motion and the rotational energy is negligible as compared to E. Hence, kj, and consequently Ue{R), do 



not depend on j in practice. Using Eq. ( 29 1 and following the developments of Appendix B, we can rewrite Eq. (|28| as 



= {2TTnf j dRdOdPdPe poiR,e,P,Pe) Pu{R,P) Pj{0,Pe)- 



(32) 



po{R,O,P,P0) and pu{R,P) are the Wigner distributions associated with ^o{R,9) and Ue{R) respectively, and 
Pj{0,Pg) has already been introduced. 

Within the harmonic approximation of Vg, po{R, 9, P, Pg) is found from Eq. (5.20) of ref. [5] to have the form 



Po{R,9,P,Pg 



(33) 
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Consequently, E'^ becomes 

^'e « (34) 

with 

h= j dRdP e-2"«(«-fl.)V?^ ^-PV(2..^n)- ^^(^^ (35) 

and 



Since Ue{R), and consequently pu(R,P), do not depend on j, the j dependence of is only due to l2- Therefore, 



Ii can be transfered into the proportionality factor of Eq. ( 34 ) , hence leading to 



^'e « (37) 



With ag — 5, corresponding to an initial angular distribution spreading over ^ 20 degrees, and 9e = 0, we arrive at 
the distribution Pj represented by the black curve in Fig. [6j This distribution is "exact" within the assumptions of 
the present Franck-Condon model. 
Besides, S-^ is equally well given by 



— 



dRde sine ^^j^*{R,e) ^t{R, 



(38) 



(see Eq. ([8|)), which is Eq. (28) at time t instead of time 0. Using Eq. (291 and following the developments of 



Appendix B, we can rewrite Eq. (38 1 as 



E^^ = i27rhy / dr ptiT) pu{R,P) Pj{0,Pe 



(39) 



where T = {R,e,P,Pe). Just as Eqs. (18) and (23) lead to Eq. (25), Eqs. (23) and (39) lead to 



E^^ = i2nhf J dTo poiTo) Pu{Rt,Pt) Pj{Ot,Pet), 



(40) 



To being the value of F = {R,9,P,Pg) at time 0. When making t tend to infinity, one recovers the semiclassical 
Wigner expression analogous to Eq. ( 25 ) in the present case where the vibration motion is frozen. 

Since the radial motion is uncoupled with the angular motion, Rt and Pt are functions of Rq, Pq and t, while 9t 
and Pet are functions of 0o, Pea and t. Using Eq. (33), Eq. (40) can thus be rewritten as 



E^^cx/i(t)/2(t) 



(41) 
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with 



I^{t) = y dRodPo e"2"H 



PuiRuPt) 



and 



(42) 



(43) 



Since pu{Rt, Pt) does not depend on j, the j dependence of S-^ is only due to hit)- Therefore, Ii{t) can be transfered 



into the proportionaUty factor of Eq. (41 1, leading thereby to 



(44) 



AB rotating freely after the photon absorption, 9t and Pg^ are given by 



= ^0 H 



P90 
mr'i 



(45) 



and 



Pet = Peo- 



(46) 



Eq. (45) is the solution of Eq. (5.4c) of ref. [6], remembering that the reduced moment of inertia / (see Eq. (A. 12)) 
reduces here to mr^, as stated at the beginning of this section. Therefore, we finally arrive at 



dOodPoQ e' 



-2ag 



Pi 



Peo 



t.Pf 



eo 



(47) 



Note that Eqs. (36) and (37) are recovered from Eq. (47 1 at time 



With the hypothetical values m = 5 g.mol ^ and re — I A, Eq. (47) applied aX t ^ 50, 100 and 500 fs leads to 



the rotational state distributions displayed in Fig. |6j in addition to the "exact" one at time 0. Beyond 500 fs, the 
distribution does not evolve. The degrading effect previously outlined is patent, the distribution getting strongly 



altered after only 50 fs. This is a clear illustration of the inability of Eq. (25), as such, to correctly describe partial 
cross sections. 



We note from Eq. ( 47 1 that the degrading effect is due to the strong variation of the rotational Wigner distribution 



along classical orbits, as previously outlined. The analogous variation being much weaker for the vibration motion, 
so is the corresponding degrading effect. 



3. Using Brown and Heller expression at a large distance rather than at a large time 



For the previous Franck-Condon process, it is equivalent to state that the distribution is exact at time 0, or at the 
time beyond which Pg ceases to vary. For a general process, the latter corresponds to the instant where the system 
crosses the frontier separating the interaction region from the free products. This frontier being well defined by a 



given value Rf oi R, we shall use instead of Brown and Heller expression (25), 
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{2Trhf J dTo po(ro) Ptr{Rf,Pf) Pn(r/,p/) Pj{ef,Pef), 



(48) 



where Pf is the value of P at i?^ when starting from Fq, with a similar definition for rf, pf, Of and Pg j. One will not 



forget to substitute 9 for 9f according to Eqs. (261 and (27) (with Of instead of Ot). The above expression is assumed 
to minimize the degrading effect, and is exact in the Franck-Condon limit. 



In practice, the delta function in ptr(Rf-, Pf) (see Eq. ( 19 )) may be replaced by a bin centered at Ti /c^ / (2^), much 



narrower than the distribution of P'j /{2^). In this work, the bin width is taken at five percent of the full-width-at- 
half-maximum (FWHM) of the previous distribution. 

As stated in the introduction, we call the present method forward I and apply it to the photodissociation of methyl 
iodide in section fVl 



E. Approximate formulation 



Eq. (48) contains a delta distribution through the translational Wigner function ptr{Rf, Pf), given by Eq. (19 1. 
This term, replaced by a narrow Gaussian or a thin box in practical calculations [^llOj. makes them heavy. However, 



at the exit of the interaction region, the classical Hamiltonian ( 24 1 reads 



He = — +E,„t{r,p,Pe) 
2fi 



with 



E^rit{r,p,Pg) 



2m 



-Ve{r) 



Pl_ 

2mr^ 



(49) 



(50) 



Eintir,p, Pe) represents the internal energy of AB. From Eqs. (15) and (49), one may thus rewrite ptr{Rf , Pf) as 



He - E + Enj - Eintirf,Pf,Pef) 



(51) 



Note the disappearence of the term 0(P/) as compared to Eq. (19), for Pf is necessarily positive. 

For not too excited rovibrational states, Enj is very well approximated by the sum of E'^, the vibrational energy 
corresponding to the n*'' state of the non rotating AB diatom, and the rotational energy Ej — h^j{j + l)/{2mrl). 
The values of ry and pf corresponding to E^ define an elliptic-like curve in the {r,p) plane (a true ellipse for a purely 
harmonic oscillator) while the values of Pgf corresponding to Ej are ~ ±?i(i + 1/2). However, rf, pf and Pef are 
weighted in Eq. (48) by pn{rf,pf) and pj{Of,Psf) which broadly extend around the previous values. Consequently, 



one expects Eint{f f ,Pf , Pe f) to have roughly the same chance to be larger or lower than E^j- In other words, their 



average difference should be negligible as compared to E, so the partial cross sections obtained by means of Eq. ( 51 ) or 



Ptr{Rf,Pf) 



2TTh 



5{He - E) 



(52) 



should not be very different. 

The interest of this approximation is that He, as a constant of motion, can be expressed in terms of Tq. One can 
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use this fact to analytically integrate with respect to Pq and po as follows. Setting 



(53) 



and 



(54) 



with [3 in the range [0,27r], Eq. (24) becomes 



He = Y + ^ + Ve{Ro,ro,eo). 



(55) 



From Eqs. (48) and (52)-(55), the partial cross section S^-* reads 



cx 



rjdrjdpdRadradBodPeo paiXo) S y ^ Q Pnirf,Pf)P]iOf,Pef) 



(56) 



with 



Q = E-^-VeiRo,ro,eo). 



(57) 



Replacing rjdri in Eq. (56) by drf /2 makes the delta function disappear and we finally arrive at the useful expression 



S'^-' cx j d/SdRodraddodPoQ /5o(ro) Pn{rf,pf) pj{df,Pgf). 



(58) 



The values of Rq, tq, 0o and Pqq contributing to the integral are those making Q positive or zero. Po and pa, which 



complete the set of initial conditions, are determined by means of Eqs. (53) and (54) with 



V = (2Q) 



1/2 



(59) 



We call the present method forward II. 

Another possibility would have been to follow Goursaud et al. jllj and integrate over one of the two momenta Pq or 
Po- But a term diverging at the boundaries of the available phase space volume would have appear in the integrand, 
rendering thereby the numerical calculation of more tricky. 
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III. BACKWARD SEMICLASSICAL WIGNER APPROACH 



From Eqs. (48) and (51), we have 



TTji = {2TThf j dFo po(ro) 5[H,-E + Er,j-E,r,t{rf.Pf.Pef)\ Pn{rf,Pf) Pj{Of,Pef) 



(60) 



In refs. [TMT5lli5H5T] . it is shown that an alternative set of coordinates to Fq is {t,He,rf,pf,6f,Pgf). The origin 
of time corresponds to the instant where the system is at Rf. The quadruplet {r f,pf,6 f, Pg j) specifies the internal 
state of AB at time and forces P/ to take the value 



Pf^[2ti{H,-EUrf,Pf,PefW/^ 



(61) 



(see Eq. (49)). Tf — {Rf,Pf,rf,pf,Of,Pgf) lies along a given trajectory. Now, any point along this trajectory can 
be reached from by moving along the trajectory a given period of time \t\ either forward (t > 0) or backward 
(t < 0). In other words, for a given Rf, {He,rf,pf,9f,Pgj:) imposes the classical path, and t the location along it. 
Consequently, {t, Hf^,rf,pf,0 f, Pg f) allows to span the whole phase space. 
In addition to that, one may show the important property p!6Hl8l I49H5T] 



dFo ~ dtdHedrfdpfdOfdPgf. 



(62) 



The exact demonstration of the above identity is not given in the previous references, but it closely follows, for 



example, the developments in Appendix C of ref. [18] for different (though partly related) coordinates. From Eqs. (60) 



and (62 1 and the straightforward integration over H^, we finally arrive at 



= (27r?i)^ / drjdpfdOfdPgf Pn{rf,Pf) Po{9},Pgf) I dt p„{Tt) 



(63) 



Integration over forces the latter to be equal to £^ — Enj + Eint{rf,Pf,Pgf)- From Eqs. (15l and (61), we thus 
have 



Pf = Tik 



nj ■ 



(64) 



At last. Ft in Eq. (63) is the value of F at time t when starting from F/ at time (the meaning of Fj is thus different 



here and in section II A ) 



To summarize, the internal state {rf,pf,6f,Pgf) of AB is randomly chosen within appropriate boundaries. Together 
with Eq. (64), they allow to generate a trajectory from Rf aX time 0. The trajectory is then propagated backward 



in time, i.e., in the direction of the reagent molecule, and po(Ft) is time-integrated until the trajectory recrosses Rf 
towards the products. The result is multiplied by the statistical weight pn{rf,pf) pj{9f,Pgf) in order to get the 



integrand of Eq. (63). A simple Monte-Carlo procedure can then be used to estimate E^-'. 



In practice, the power of the backward approach is limited by the fact that one cannot a priori guess which values 
of rf, Pf, Of and Pgf lead to trajectories crossing the Wigner region, corresponding to the phase space volume 
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where Wigner distribution po{T) takes significant values (there is some arbitrariness in this definition). Therefore, 
a straightforward application of Eq. (63) may require running a large amount of useless trajectories that do not 
contribute to S^^* . To go round this difficulty, one may first apply the forward approach in order to determine the 
boundaries of rj, pf, Of and Pg f contributing to S^"', and then apply the backward method with these variables 
selected within the previous boundaries. 

The practical method used here to perform this selection is as follows. First, one runs a few thousand trajectories, 
say N, within the forward I method, leading to the same number of final points (r^, p/j, 9k, Psk): ^ = l,...,iV, at 
R — Rf. Next, one randomly generates a point (ry, pf, Of, Pe f) and checks whether it lies within at least one of the 
small rectangular cuboids defined by |r/ - rfc| < 77^, \pf - Pk\ < Vp, 1^/ - ^'fel < Ve, \Pe / - -Pe^l < ?7P(,, fc = 1, ^• 
If so, this point serves as initial conditions together with R = Rf and Eq. (64). Otherwise, one randomly generates 
another point and so on. The parameters rjr, ^jp, rjg and ijpg have to be chosen from a visual inspection of the domain 
covered by the rectangular projections of the cuboids in the planes (r,p) and {6, Pg), separately. The parameters must 
be large enough for the domains to be compact, as they appear to be when running millions of trajectories within the 
forward I method. 

Finally, we have found that for the model of methyl iodide photodissociation considered in section [V] the values 
of knj appear to be almost independent on j (see section |II D| for an explanation), a bit less on n. In such a case, a 



single batch of trajectory can be run with Pf — k„o for calculating all the s corresponding to n 



IV. THE STANDARD METHOD OF GOURSAUD ET AL. 

The method of Goursaud et al. was initially applied to a bi-dimensional model of triatomic ion fragmentation 
with frozen valence angle. The method was later extended by Schinke to realistic three-dimensional triatomic 
photodissociations [48j. Within this approach, which is detailed in Chapter 5 of Schinke's book [B], E^-* is given by 

^E^i^ + Sjo) J dTopo(To)sineo6{He-E)6{nirf,pf)-n)di\Pef\-j) (65) 
where the only quantity not defined until now is the final vibrational action n(rf,pf) of AB, given by 



Arf,Pf) 



dr p[r) 



(66) 



with 



p{r) = (2m E,ntirf,pf,Pef) 



1/2 



(67) 



and Tout are the values of r at the inner and outer turning points. The internal energy Emtir f , p f , Pe f) is given 
by Eq. 

In practice, the delta functions are replaced by standard bins (SB) of unit height and width, or Gaussians the 
FWHM of which is usually taken at 10%. The second procedure is called Gaussian binning (GB) [S2HSH]- These 
procedures are discussed at length in ref. [T5]. Since no Wigner distributions are used to weight the final dynamical 
states, the present method is more classical than the previous ones. 

The transition moment being absorbed in pq{Tq), Eq. (65) is similar to Eq. (5.23) in ref. [6] (see also Eq. (5.22)). 
The main difference is that the degeneracy factor (1 + Sjo) has been added in the present work. Note that Schinke 
and co-workers also include this factor in practice jSUI . 

The degeneracy factor doubles the integral in Eq. (65) for j = 0. This counterbalances the fact that due to the 
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(5(|Pe/| — j) term, two values of Pgf contribute to the integral for j > (±j), against only one for j = 0. The interest 
of this factor clearly appears in the purely statistical limit where n(rf,pf) and Pgf are random variables. S^-* being 
proportional to the density of probability to get n{rf,pf) = n and \Pef \ ~ j (see Eq. ([65|), all the S^'^'s are equal 
(without the degeneracy factor, E^*^ would be half Yi^'^^), a result in conformity with the quantum phase space 
theory expectations of equal final state populations |BDJ [HI] • 

V. PHOTODISSOCIATION OF METHYL IODIDE 

We briefly summarize the main features of the reduced-dimensionality model of methyl iodide photodissociation 
and the quantum dynamical method used to check the validity of the semiclassical Wigner method. More details can 
be found in the references mentioned below. 

A. Model 

The CH3I molecule is considered as a CXI pseudotriatomic molecule [1TJ|S11[H3], the pseudoatom X=H3 being 
located at the center-of-mass (CM) of the three H atoms. R is the vector between the CH3 (or C— X) CM and I and 
r is the one between X and C. r represents the umbrella bend of the C— H3 group. 

Photodissociation of CH3I is assumed to take place upon optical excitation at 266 nm (A band) from the X^Ai 
ground electronic state to the ^Qo, ^Qi and ^Qi excited electronic states. Taking however into account that absorption 
to the ^Qi state is relatively small at 266 nm, the present simulations only involve the ^Qo, and ^Qi excited electronic 
states, in addition to the ground state X^Ai. The X^Ai and ^Qi electronic surfaces correlate asymptotically with 
the CH3 + I(^P3/2) products, while the '^Qo surface correlates with the CH3 + I*(^Pi/2) products. In addition to the 
coupling of X^Ai to the excited electronic states through electric-dipole moments, the '^Qo and ^Qi states are non 
adiabatically coupled. Taking into account the transitions between these states in the semiclassical Wigner method 
would require treating these by means of a semiclassical approach of non adiabatic transitions such as, for instance, 
the Landau-Zener model |M| , the Zhu-Nakamura model [55] , or the surface hopping method of Tully [551 E2 ■ For 
clarity's sake, however, we artificially take at zero the coupling between the '^Qo and ^Qi states, focusing our attention 
on the semiclassical Wigner method for fragmentations taking place on a single excited electronic state. However, we 
plan to extend this method to processes involving non adiabatic transitions in a near future. 

Upon optical excitation at 266 nm, the energies available to the final products in the '^Qo and ^Qi states are 
11258.53 cm-1 and 18862.09 cm-1, respectively 

High-quality ab initio calculations have been used to model the three electronic potential energy surfaces (PESs) 
involved in the calculations. In the case of the X^Ai ground state, the PES {Vg{R,r,9)) is represented as a sum of 
three potential interactions in the Rc-i (the C— I nuclear distance), r and 9 coordinates, respectively. The interaction 
potential in the Rc-i coordinate is taken from the recently reported 2D ground-state potential for CH3I, obtained by 
means of multireference spin-orbit configuration interaction ab initio calculations |68l I69j . The potential interactions 
in the r and 9 coordinates are represented by harmonic oscillator functions [55J HI] . 

The PESs for the excited electronic states '^Qq and ^Qi {Ve{R,r,9)), are the ab initio PES constructed by Xie et 
al. [70]. These are improved versions of the previous nine-dimensional surfaces of Amatatsu et al. [71], where the 
spin-orbit configuration interaction method was used with a better basis by changing the valence double-^ level to the 
triple-C one. Out of the nine coordinates of the surfaces, the six coordinates which are neglected in the simulations 
were fixed at their equilibrium values |62j . 

The electric-dipole moment functions {d{R,r,9)) coupling radiatively the X^Ai ground state with ^Qo and ^Qi 
were obtained from the 116 initio calculations of Alekseyev et al. [521 IZl] , subsequently fitted to analytical forms used 
in the simulations [62j. 
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B. Time-dependent quantum calculations 



CH3I is initially in the rovibrational ground state (/)o(i?, r, 0), variationally obtained from Eq. (131 within the 



framework of an adiabatic approximation detailed in ref. |26j . CH3I is then excited to one of the '^Qo a-nd ^Qi 



electronic states, creating a wave packet ^o{R,i^,9) (see Eq. (12)) undergoing dynamical evolution according to 



Eq. (10). In order to solve this equation, the wave packet $t(i?, r, 6*) is represented in a basis set consisting of a 
two-dimensional rectangular grid for the radial coordinates and an angular basis including 24 Legendre polynomials 
for the 9 coordinate. The rectangular grid consists of 450 equidistant points in the R coordinate in the range 
3.5 flo < -R < 16.0 flQ, and 32 equidistant points in the r coordinate in the range — 1.6ao <r< 1.6 oq. Propagation of 
the wave packet is performed by representing the evolution operator by means of a Chebychev polynomial expansion. 
The wave packet was propagated for 200 fs with a time step At = 0.4 fs, and was absorbed at the edge of the grid 
in the R coordinate after each propagation time step by multiplying it by the function exp[—a{R — Rabs)^], with 
a — 0.9 Oq ^ and Rabs = 13.0 qq. In order to obtain the product fragment distributions of interest, the wave packet 
was projected out in the asymptotic region onto the fragment states. Details of the projection procedure are given in 
ref. [26]. 



C. Comparison between semiclassical and quantum results 

$0(^1 fi Q) was expressed as a product of three independent Gaussians respectively depending on R, r and 9, a very 



good approximation in the present case. As in section II D the density po(ro) was then deduced from Eq. (5.20) of 
ref. 0. 

The results obtained by means of the forward I, forward II and Goursaud et al. methods involved 5 million, 100 
thousand and 1 million trajectories, respectively. The results obtained by means of the backward method involved 60 



thousand trajectories per vibrational level (see end of section III). These numbers are for each electronic state, Qo 
and ^Qi- For the backward approach, 77^ and rjg were kept at 0.02 while rjp and rjpg were taken at 0.2 (in atomic units 
for rjr. Tip and ryp^ and in radian for rio). Rf was found to be equal to 13 and 10 bohr for the ^Qq and ^Qi states, 
respectively. 

The vibrational state populations P„, deduced from Eq. (|3| by summing over j, are given in Fig.[7] The vibrationally 
resolved rotational state distributions are displayed in Fig. [8] and Fig. [9] The latter are simply denoted Pj , but the 
value of n to which they refer is indicated. 

The agreement between backward and quantum mechanical (QM) results is quantitative for all the distributions 
but Pn in the ^Qi state, for which it is nevertheless very good. 

These conclusions hold between forward I and QM results, apart from Pj, n = 1, 2, in the ^Qi state for which the 
quality of the agreement decreases. The reason seems to be the following. As seen at the end of section |IID[ the 



partial cross section involves a delta function through ptr{Rf,Pf) (see Eqs. (48 1 and ( [T9| ) which is replaced by a 
narrow bin. However, the bin cannot be too narrow, for a negligible amount of trajectories would contribute to the 
partial cross sections. With a true delta function, P would be strictly equal to hknj- With a bin, however, there is 
an uncertainty on the value of P around hk„j, which seems to cause the differences observed between the forward 
I results and the backward and quantum ones. By dividing the width of the bin by two, we indeed reduced the 
differences. We tried to reduce even more the width of the bin, but this prevented from converging the partial cross 
sections. The minimum number of trajectories necessary to converge the calculations was found to be 800 thousand 
for the forward I method in both electronic states, against 10 thousand for n = in the '^Qq state, 20 thousand for 
n = 1 in the ^Qq state, 8 thousand for n = in the ^Qi state, 10 thousand for n — 1 in the ^Qi state and 30 thousand 
for 71 = 2 in the ^Qi state. As a matter of fact, the backward method is at the same time more efficient and more 
accurate than the forward I method. 

The agreement between forward II and QM results is surprisingly good for all the rotational distributions (we 
have no explanation for that), a bit less for the vibrational ones, especially in the ^Qi state. Like the forward I and 
backward methods, the forward II method accounts for the vibrational inversion in the ^Qi state. 
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The agreement between Goursaud et al. and QM results is very satisfying for the rotational distributions in the 
^Qi state, except for j — 0, and for the vibrational distribution in the ^Qq state. On the other hand, the method 
of Goursaud et al. fails at reproducing the vibrational inversion in the ^Qi state, and badly describes the rotational 
distributions in the '^Qo state, even qualitatively. Overall, this method is only semi-quantitative. 

As previously stated, the energy available to the final products is larger in the ^Qi than in the '^Qo state. This 
explains in part why there is more vibrational and rotational excitation in ^Qi than in '^Qo- The system is thus 
less quantum-like in ^Qi than in '^Qoj justifying thereby why Goursaud et al. method, the more classical of the four 
semiclassical approaches considered in this work, is able to reproduce the energy partitioning in ^Qi and not in '^Qq. 

VI. CONCLUSION 

In their concluding remarks. Brown and Heller W raised the basic issue of including rotations in their semiclassical 
Wigner description of photodissociation dynamics so as to make it full-dimcnsional and thus, realistic. This issue has 
been solved in the present work. 

Three methods have been proposed, respectively called forward I, forward II and backward. While forward ap- 
proaches involve trajectories run from the reagent molecule onto the products, the backward one deals with trajectories 
run in the reverse direction. This makes the backward method state-selective, and thus, ideally suited to the descrip- 
tion of state-correlated distributions measured in velocity imaging experiments. 

The forward I and backward methods are exact applications of the semiclassical Wigner treatment. They closely 
parallel the quantum description as far as possible. The forward II approach involves an approximation strongly 
increasing the numerical efficiency of the semiclassical Wigner treatment as compared to the forward I method. 

In addition to the usual vibrational Wigner distributions, these three approaches include rotational Wigner dis- 
tributions which seem to be introduced for the first time in the present context. These densities have remarkable 
structures clearly showing that classical trajectories most contributing to rotational state j are those reaching the 
products with a rotational angular momentum close to -\- 1)]^^^ (in Ti unit). 

The results obtained by means of these methods are compared with rigorous quantum results in the case of Guo's 
triatomic-like model of methyl iodide photodissociation |41) and very good - often quantitative - agreement is found, 
especially with the forward I and backward methods. In comparison, the standard and more classical method of 
Goursaud et al. |11) is only semi-quantitative. Last but not least, the backward approach appears to be much more 
powerful, and even more accurate than the forward I method which requires far more trajectories than the former 
approach to provide converged results. 

This study demonstrates the applicability of the semiclassical Wigner treatment to realistic triatomic photodisso- 
ciations and confirms its level of accuracy as compared to the initial work of Brown and Heller [3] . Important next 
steps are the extention of the method to processes involving non adiabatic transitions and/or polyatomic species. 
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Appendix A: Derivation of Eq. (|8| 

Consider a given state ^'•■'*^'(R, r) of ABC for the value J of the total angular momentum quantum number, the 
value M of its projection on the Z-axis of the laboratory reference frame, and the parity e under inversion of R and 
r. This state can generally be expanded as [531 121] 



(A.l) 



with 



2J+ 1 



167r2(l + ,5f,o) 



1/2 



[Dii.ni^, 7) + e(-l)('^+'')Z?i. /?, 7)]- 



(A.2) 



a, /3 and 7 are the Euler angles orienting ABC in the laboratory reference frame [71]. The helicity quantum number 
n is the projection of the total angular momentum on R, chosen as the Z-axis of the body-fixed frame. Dj^j ^(a, /3, 7) 
and -D}{^ _f2(a, /?, 7) are Wigner D-matrix elements [71]. 

When J — Q, M = Vl = {) and e is necessarily equal to 1 for VF//fj(Q;, /3, 7) no to be 0. Given that q(q;, /3, 7) = 1 



(see Eq. (4.1.26) of 74 ), Eq. (A.l) reduces to 



(A.3) 



We can thus rewrite \I'^^ (R, r) and $t(R, r) in Eq. as 



(A.4) 



and 



$,(R r) = ^^^^''^'^^ 



(A.5) 



where for convenience's sake, subscripts and indices relative to angular momentum quantum numbers and parity have 
been dropped. 

Eqs. 0, ( |A.4[ ) and ( |A.5| ) lead to 



(A.6) 



Moreover, some steps of algebra allow to prove the two identities 



dRdr — R r dRdrsinOdOdasinfidpd'y 



(A.7) 
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and 



/ 



dasinfidpd"f = Stt^ 



(Ai 



Eqs. (|X6t-(|X8l) finally lead to Eq. 

For J = 0, ^'^■'(R, r) and $t(R, r) are shown to satisfy the Schrodinger equations [75l [76] 



(A.9) 



and 



at 



(A.10) 



with 



T 



d „o d d ^ d d . B 



smb 



2/xi?2 ai? dR 2mr^ dr dr 2Isin9 86 86 



(A.ll) 



The reduced moment of inertia / satisfies the identity 



1 _ 1 1 



(A.12) 



From Eqs. ( [XI] ), ( [Xsj ), (|X9|-(|ATT]), ■i''^\R,r,e) and $t(i?,r,6l) are shown to satisfy Eqs. (|9|-([TT]) 
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Appendix B: Derivation of Eq. (18 1 

Any overlap S defined by 

E = y dRdrd0^l{R,r,0)^2{R,r,( 
can be rewritten as [U [TOl [12] 

S = (27r?i)3 y dRdrdOdPdpdPg pi{R, r, 9, P,p, Pg) p2{R, r, 9, P,p, Pg) 



(B.l) 



(B.2) 



where pi{R, r, 9, P,p, Pg), I — 1 or 2, is the Wigner density defined by 

p,{R,r,9,P,p,Pe) = j-^ J dSRdSrdsg e'KPsH+pSr+Pese)/h + + + ^^(^R^SR,r-Sr,9-Sg). (B.3) 

This expression is a generahzation of Eq. ([T]) to three dimensions. 



A pedestrian demonstration of the strict equivalence between Eqs. (B.l) and (B.2) for one configuration space 
coordinate is given in Appendix B of ref. |18j . In the present case of three coordinates, the developments are more 
tedious, but present no difficulty. 

Setting 



*i(i?,r,( 



71J 



1/2 



and 



^2{R-r,9) = ^t{R,r,9), 



(B.4) 



(B.5) 



we arrive from Eqs. (B.2) and (B.3I at 



YTJ = {2T:nf lim / dRdrd9dPdpdPe pt{R, r, 9, P,p, Pg) ptr{R, P) Pn{r,p) Pj{9, Pg). (B.6) 



Pt{R, r, 9, P,p, Pq) is related to ^t{R, 0) by Eq. ( |B.3[ ). The translational Wigner distribution ptr(R, P) is given by 

(B.7) 



Ptr{R, P) 



j dse^'P'"' 


\ ^ 1 







while the vibrational and rotational Wigner distributions Pnir,p) and pj{9,Pg) are given by Eqs. (20) and (21 1 



Eq. (B.7 1 gives 
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which, using 



S{x) 



2tt 



(B.9) 



leads to 



PtriR,P) 



S[2{p-hk„j)/h] 



(B.IO) 



or equivalently, Eq. (19). Eq. (B.IO) is indeed readily obtained from Eq. (19) by means of the following theorem 



(B.ll) 



where the Xk^s are solutions of f{x) = [75] . 
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Appendix C: Semiclassical limit of the rotational Wigner distribution 



Our goal here is to derive an analytical expression of pj{6,Pg) within the semiclassical approximation along the 
relevant directions defined by = 7r/2 and Pg ~ 0, ±{j + 1/2). Pg is expressed in h unit in throughout this appendix. 



These directions are emphasized in Fig. 10 for j = 10. The first two ones (Di and D2) are the two orthogonal 
symmetry axis of Pj{0, Pg). The next two ones (D3 and D4) are not exact symmetry axis, but they can be considered 
as local symmetry axis of the lateral ridges for sufficiently large j's, as shown further below. 
The semiclassical (WKB) limit of the Legendre polynomial Pj{cos9) is given by [77] 



Pj{cos6 



2cos[{j + 1/2)9 -tt/A] 
[27r(j + l/2)sm6l]i/2 ' 



(C.l) 



From Eqs. (16 1, (21) and (C.l), we have 



Pj{e,Pg) = ^ j dse^'^o' [sin{e + s)sin{e - s)Y''^ cos [(j + 1/2) (61 + s) - 7r/4] cos[{j + l/'2){e - s) - t: /A]. (C.2) 



Using the fact that 



cosx 



(C.3) 



we arrive after some steps of simple algebra at 



,Pe) = ^ / dse^'^o" [sin{9 + s)sin{e - s)Y'^ [cos[{2j + l)s\ + sin[{2j + 1)9]] 



(C.4) 



Along direction Di, defined by ^ = 7r/2, we have 



Pj{^l2,Pe) = ^ ds e^'Po' coss [cos[{2j + l)s] + (-1)^"]. 

'^'^ J—k/2 



(C.5) 



The boundaries of the integral are determined by Eq. (22). From Eq. (C.3) and given that 



7r/2 



ds 



siniira) 



7r/2 



Trsinc{Tra), 



(C.6) 



we finally obtain 



Pj{n/2,Pg) = ^(-IF [sinc[n{Pg + 1/2)] + sinc[Tr{Pg - 1/2)] 

+ [sinc[TT{Pg + j + 1)] + sinc[7r{Pg + j)] 
+ sinc[Tr{Pg — j)] + sinc[TT{Pg — j ~ 1)]] . 



(C.7) 
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This expression is represented in Fig. [TT] for j = 0,1,2 and 6 and is compared with the exact Wigner distribution. 
For j = 0, the agreement is aheady satisfying. For j — 1 and 2, it is very good and for j = 6, there is virtuaUy no 
difference between both results. 



For large j's, the central peak is due to the first pair of sine functions in Eq. (C.7) and its sign is thus given by 
(— l)-'. The lateral peaks, or cuts of ridges, centered at —{j + 1/2) and (j + 1/2) are due to the second and third 
pairs of sine functions, respectively. They are always positive. Each of these three peaks is thus made of two narrower 
peaks sufRciently close to each other for resulting in a single peak. The full width at half maximum is 1.2 for the sine 
function and 1.64 for the central and lateral peaks. 

For small j's, the situation is a bit more complex, since the three peaks overlap. Consequently, the lateral peaks 
are not exactly centered at —{j + 1/2) and (j + 1/2). For j = 0, the lateral peaks have merged with the central peak, 
and their top is rigorously at 0. 



Along direction D2, defined by Pg — 0, we have from Eq. (C.4) 



Pj{e,0) = ^ / ds [sin{e + s)sin{9 - s)]^/^ [cos[{2j + l)s] + sin[{2j + 1)6]] (C.8) 



where the boundaries are found from Eq. (22) to be given by s± = ±6 if < 7r/2 and s± = ±{t: — 9) ii 9 > tt /2. It is 
thus clear that the angular range around 7r/2 mostly contributes to pj{9,0). 

For significantly large j's, cos[{2j + l)s] strongly varies within the previous range and can thus be neglected with 
respect to sin[{2j + 1)9]. We thus arrive at 



p.[0^O) ^ —sin[{2j + 1)0]] I ds[sin{9 + s)sin{9- s)]^/^. (C.9) 

The integral seems not to be analytically calculable, but it appears that a good approximation of it is 2sin9^^^. 
Therefore, pj{9,0) reads 



{9, 0) w ^sin9^^^sin[{2j + 1)9]] . (C.IO) 



This expression is represented in Fig. [12] for j = 0,1,2 and 6 and is compared with the exact Wigner distribution. 
The agreement is correct for j = and 1, good for j = 2 and very good for j = 6. 



Along direction D3, defined by Po = j + 1/2 for sufficiently large j's, we obtain from Eqs. (C.3| and (C.4 1 



Pj{9J + 1/2) [ ds [sin{9 + s)sm{9 - s)]^/^ 



\ I ds e2^(2j+i). + _ s)]i/2 (c.ll) 



47^3 



-\ [ ds e*(2j+i)s rgjj^(-g) ^_ s)sin{9 - s)]^/'^ sin\{2j + 1)9] . 
27r3 7 



Since the second and third integrands strongly oscillates as compared to the first one, the above expression is well 
approximated by its first line. Using the previous approximation for the integral of the first line, we then arrive at 
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P.-(^?,J + l/2) 



95/2 



27r3 



(C.12) 



The quality of the agreement between this expression and the exact Wigner distribution is shown in Fig. [13] for j = 10. 

Finahy, the semiclassical expression of the rotational Wigner distribution along direction D4, defined by Pg = 
— {j + 1/2), is also given by the right-hand-side of Eq. (C.12). due to the symmetry of the density with respect to D2. 
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FIG. 1: Blue curve: r-dependence, in the ^Qi electronic state, of the potential energy in the separated products {R > 12 
bohr) for the reduced dimensionality model of methyl iodide photodissociation of Guo [41]. Its shape is the same in the ^Qo 
electronic state, with the difference that the bottom of the well is 7603.92 cm~^ higher. Red curve: second order development 
of the previous potential. 





FIG. 2: Perspective views of the vibrational Wigner distributions for the levels n = — 5 for the potential energy shown in 
Fig. [l] (blue curve). The r and p axis are directed towards the right and left, respectively, r belongs to the range [-1.2,1.2], and 
Pr to the range [-12,12], in atomic units. 




FIG. 3; Perspective views of the rotational Wigner distributions for the states j = — 3,6 and 10. The 9 and Pq axis are 
oriented towards the right and left, respectively. 6 belongs to the range [0,7r], and Pe to the range [-15?i,15?i]. 




FIG. 4: Front view of p(i,{9,Pe) along the 6'-axis, directed towards us. The wells along the symmetry axis of the distribution 
correspond to negative peaks, not visible in Fig.jsj having magnitudes comparable with those of the positive peaks. They form 
the black downward peak, resulting from the alignment of 6 peaks. The two lateral peaks, forming ridges along the S-axis, have 
been found to be approximately centered at ±[j(j + 1)]^'^^ in Tt- unit. 
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j 



FIG. 5: The difference between (i) the exact value of Pe corresponding to the summit of the rotational Wigner density in the 

upper half plane and (ii) the quantum value h[j{j + 1)]^^^, tends to as j takes large values. In the particular case where 
j = 0, the lateral ridges have merged with the central peak, and their top is rigorously at the expected quantum value, i.e., 0. 




FIG. 6: 



Time dependence of the rotational state distribution for the Franck-Condon process. Time is given in fs. 
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FIG. 7: Vibrational state distributions in the ^Qo and ^Qi states found by means of the QM, forward I, forward II, backward 
and Goursaud et al. methods. 




FIG. 9: Vibrationally resolved rotational state distributions in the Qi state found by means of the QM, forward I, forward II, 
backward and Goursaud et al. methods. 
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FIG. 10: Contour plot representation of pio{0,Pe) and directions Di, i = 1 — 4, along which semiclassical expressions of the 
rotational Wigner distribution are derived in Appendix C. Po is in h unit. Solid and dashed contours correspond to positive 
and negative densities, respectively. 




FIG. 11: Comparison between the semiclassical value of pj{-K /2, Pe) along the Di direction (see Fig. |10[ ), obtained from 
Eq. (C.7l, and its exact quantum value for j = 0, 1, 2 and 6. Pe is in Ti unit. 




FIG. 12: Comparison between the semiclassical value of pj{9,0) along the D2 direction (see Fig. |10[ ), obtained from Eq. ( |C.10[ |, 
and its exact quantum value for j = 0, 1, 2 and 6. 




FIG. 13: Comparison between the semiclassical value of pj{6,j + 1/2) along the D3 direction (see Fig. |10[ ), obtained from 
Eq. (C.12I, and its exact quantum value for j = 10. 



